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Abstract. In this paper we construct the quantum spectral curve for the quantum dynami- 
cal elliptic 0[„ Gaudin model. We realize it considering a commutative family corresponding 
to the Felder's elliptic quantum group i?T-,B(0^ri) ^.nd taking the appropriate limit. The 
approach of Manin matrices here suits well to the problem of constructing the generation 
function of commuting elements which plays an important role in SoV and Langlands con- 
cept. 

Key words: Manin matrices; L-operators; elliptic Felder i?-matrix; Gaudin models 
2000 Mathematics Subject Classification: 37K15 



1 Introduction 

The Gaudin model plays intriguingly important role in the modern mathematical physics as 
like as in purely mathematical subjects as the Geometric Langlands correspondence. In fact, it 
is shown that the separation of variables of the quantum rational Gaudin model is equivalent 
to the categorical part of the geometric Langlands correspondence over the rational curve with 
punctures over C. Its physical importance lies in the condense matter field, this system pro- 
vides an example of the spin interacting magnetic chain. This paper deals with a new useful 
formalism applicable to the dynamical elliptic case of the Gaudin model. This case provides an 
interpretation for the Langlands correspondence over an elliptic curve. From the physical point 
of view this case is responsible for the periodic spin chains. 

The main strategy of the paper goes throw the classical ideas of [H [2] . We start with the 
formalism of L-operators corresponding to the Felder "elliptic quantum groups" ET-fi{Q\^) [Sll^. 
Using the iiLL-relations we construct the Bethe elliptic commutative subalgebra from these 
L-operators. To make the proof of commutativity more apparent we use technique of Manin 
matrices developed in [5l E]. For some particular cases this commutative subalgebra can be 
found in [7] in a slightly different form. It is worth to mention the detailed work [8] described 
the centre of ET-fi{5{n)- We refer also to the rational version in [9l [10], to the trigonometric 

*This paper is a contribution to the Proceedings of the Workshop "Elliptic Integrable Systems, Isomonodromy 
Problems, and Hypergeometric Functions" (July 21-25, 2008, MPIM, Bonn, Germany). The full collection is 
available at |http: / /www. emis.de/journals / SIGMA /EUiptic-Integrable-Systems.html , 
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dynamical case in [T2] . Then we degenerate these families to obtain a commutative family 
for the elliptic gl„ Gaudin mode0- To do it we consider a quantum characteristic polynomial, 
which has a suitable form for the degeneration. This approach extends the method applied in 
the rational case in |13| . 

The article is organized as follows: in Subsection 12.11 we consider the dynamical elliptic L- 
operators defined by the RLL relation with the Felder elliptic i?-matrix. Then in Subsection l2.2l 
we define Manin matrices following [5l [6] and demonstrate that the L-operator multiplied by 
some shift operator is a Manin matrix. In Subsection 12.31 we consider a fusion procedure of the 
L-operators and construct a commutative family as traces of the "fused" L-operators. In Subsec- 
tion [23] we show that a characteristic polynomial of the considered Manin matrix generates the 
obtained commutative family. Subsection 12.51 is devoted to commutative families corresponding 
to the traces of powers of Manin matrices related to the families constructed in 12.31 via Newton 
identities. In Subsection 12.61 we consider briefly a trigonometric degeneration of the characte- 
ristic polynomial and of the commutative families. We discuss their connection with the Hopf 
algebra f/g(0l„). 

In Subsection 13. 1 1 we obtain a characteristic polynomial which generates a commutative family 
by passing to the limit h ^ 0. In Subsection 13 . 21 we consider L-operators that gives commutative 
families for the elliptic 0l„ Gaudin model. Subsection 13.31 is devoted to a twist relating these 
L-operators with standard L-operators for the elliptic 0(„ Gaudin model. We consider 5(2 case 
as an example in Subsection 13.41 In Subsection 13.51 we present the result of application of the 
Newton identities to the elliptic gl^ Gaudin model. 

2 Elliptic quantum groups and commutative families 

Here we use the g[„ dynamical RLL relation to construct the commutative families for dynamical 
L-operators. We consider an arbitrary L-operator that defines representation of the "elliptic 
quantum group" Er^hid^n) introduced in [U 3]. The commutativity in the dynamical case is 
understood modulo Cartan elements. Restricting these "commutative" families to the space 
annihilating by Cartan elements one obtains a corresponding integrable system. 

We start with a definition of the odd Riemann theta-function related with an elliptic curve 
defining the corresponding i?-matrix. Let r G C, Imr > be a module of the elliptic curve C/F, 
where T = Z + tZ is a period lattice. The odd theta function 6{u) = —9{—u) is defined as 
a holomorphic function on C with the properties 

e{u + 1) = -e{u), 9{u + r) = -e-2™-^*^^(n), 9'{0) = 1. 

We use the matrix (or "Leningrad") notations. Let T = tj ■ aij (g) • • • (8) ajyj be a tensor 
over a ring 9\ (a complex algebra or the field C), where tj G and aij belong to a space End C". 
Then the tensor T^'^1'-'''n) the following element of 9^ (End C"j®^ for some M ^ N: 

rp{ku...,kM) = ^tj aij • • • (g) aN,j • • • (g) 1, 

j 

where each element Uij is placed to the /cj-th tensor factor, the numbers ki are pairwise different 
and 1 ^ fcj ^ M. (The condition ki < ■ ■ ■ < kj\f is not implied.) 

We also use the following notation. Let L(A) = -F(Ai, . . . , A„) be a function of n parame- 
ters Xk taking values in an algebra 21: that is F : C"' — s- 21. Then we define 

F{X + P)=F{Xi+P^,...,Xn + Pn) 



^Sometimes we omit the word "dynamical" for briefness. 
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for some P = (Pi, . . . ,i-*n), -ffc S 2t. We omit here the convergence question considering only 
such situations where this is the case. 

2.1 Elliptic dynamic i?LX-relations 

First of all we introduce a notion of a dynamical elliptic L-operator corresponding to the Felder 
i?-matrix. Using a lemma about the products of these L-operators and the particular choice of 
the L-operator we prove the commutativity of the family of operators under consideration. 

Let {ej} be a standard basis of C" and {Eij} be a standard basis of EndC", that is 
EijCk = (5^ej. In [31 H] Felder introduce the following element of EndC" ® EndC" depending 
meromorphically on the spectral parameter u and n dynamic parameters Ai, . . . , A„: 

e{u + h) 



e{u) ^ 

'Oih^ e{u - \,)e{h ) 



where Xij = Aj — Aj. It is called the Felder's dynamical R-matrix. It satisfies the dynamical 
Yang-Baxter equation (DYBE) 

R^^'^\ui-U2]X)R^^^\ui-u-i;X + hE^'^'^)R^'^^\u2-u^;X) 

= R^^''^\u2-us;X + hE^^'^)R^^^\ui-us;X)R^^^\ui-U2;X + hE^^^) (2.3) 

and the relations 

i?(-)(-n;A)i?(-)(.;A) = ^(^i±|^^, (2.4) 

(E« + ES^)R{n; A) = R{u; A) (4^^ + E^^) , (2.5) 
(51^) + 5f A) = i?(^; A) (5 W + 5f ) , (2.6) 

where Dx = J2 ^kkm^, D\ = ^ ^klwT- should comment that we always mean by 

k=i k=i * 

A the vector Ai, . . . , A„ and the expression of the type A + hE^'^^ as the argument of (12. Ih with 

R = hE\- \ The relation (1231) is obvious. The relation (ESI) follows from (|23D. 



Let 9^ be a C[[/i]]-algebra, L(n; A) be an invertible n x n matrix over $H depending on the 
spectral parameter u and n dynamical parameters Ai, . . . , A^. Let hi, . . . ,hn be a set of some 
pairwise commuting elements of . If the matrix L(n; A) satisfies dynamical RLL-relation 

_ y. X)lW (n; A + )^(2) ^) 

= L(2) (t;; A + hE^^'^)L^^'^ {u; X)R^^^^ {u - v; X + hh), (2.7) 
{Eu + /ii)L('u; A) = L(n; A)(£;ii + /i^), (2.8) 

then it is called a dynamical elliptic L-operator with Cartan elements h^. The argument X + hh 
is always meant in the sense of (12. ip with Pj = hhi. 

Let us introduce an equivalent but more symmetric form of RLL relations. For each L- 
operator we introduce the following operator (similar to an operator introduced in |llj): 

LD{u) = e~^^^L{u;X)- (2.9) 
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The equation (j2.7p can be rewritten in terms of this operator: 

R^^'^\u-v-X)L^^\u)L%\v)=L%\ij)L'^^\u)R'^'^'^\u-v-\ + hh). (2.10) 

Lemma 1. If Li{u; X) € End(C")®9^i and L2{u; X) € End(C") (8)9^2 are two dynamical elliptic 
L-operators subject to the two sets of Cartan elements: = {h\, . . . , h\) and h'^ = (/if, . . . , hf^) 
then their matrix product L2{u; A)Li(n; X+hh'^) S End(C"')(8)9^i<8)9^2 is also a dynamical elliptic 
L-operator with Cartan elements h = + = {h\ + h\, . . . , /i^ + /i^). Thus, if Li(u; A), . . . , 
Lm{u; X) are dynamical elliptic L-operators with Cartan elements h^, . . . , /i™ then the matrix 

^ f m ^ 

Yl Lj iu;X + h ^ 

m^j^l \ l=j+l 

is a dynamical elliptic L-operator with Cartan elements h = ^ /i*. 

1=1 

Remark 1. The arrow in the product means the order of the factors with respect to the growing 
index value: the expression Y\ means A2,A2Ai. 

The basic example of the dynamical elliptic L-operator is the dynamical Felder i?-matrix: 
L{u) = R{u — v;X). In this case the second space End(C"') plays the role of the algebra fH. 
Here f is a fixed complex number and the Cartan elements are /i^ = i?^^ . Lemma [J allows to 
generalize this example: let vi,...,Vm be fixed numbers, then the matrix 

^ / m 

mS^\u-{vj]-X)= J] R^^^^ \u-Vj-X + hY, E^^^ 

m^j^l \ l=j+l 

(I) 

is a dynamical elliptic L-operator with Cartan elements = ^kk- 

1=1 

The relation (12. 4p gives another important example of the dynamical elliptic L-operator - the 
matrix L(n) = R^'^^^v — u)~^ with the second space EndC" considering as $H (the number v is 
fixed), or equivalently, the matrix L{v) = R^^'^\u — v)~^ with the first space EndC" considering 
as 9^. Let us introduce the notation 

M(°)({nJ;7;;A) = R^''Uui-v;X + hf2E^'A, 

l^i^m \ l=i+l / 

where ui, . . . ,Um are fixed numbers. Then the matrix M.^^\{ui};v; X)~^ is a dynamical elliptic 

(I) 

L-operator with the spectral parameter v and the Cartan elements = ^kk- 

1=1 

A more general class of the dynamical elliptic L-operators associated with small elliptic 
quantum group erfi{Q\-n) was constructed in the work [T^. This is an C[[/i]]((Ai, . . . , A„))-algebra 
generated by the elements tij and satisfying the commutation relations 

tijhk = {hk — Sik + Sjk)tij, 
h~^(tijXk — (Afc — h5ik)tij) = 0, 

^ (j'ij^ik ^ik^ij) — Oj 

h tiktjk 7TT tjktik — 0, i T 
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where t. 



6ij+hiij, xjj^ 



,{2} 



0, iT^k,j^ I, 



Xi — Xj , X]j = Xi — Xj — hhi + hhj and the elements hi, . . . ,hk, Xi, . . ., 
Afc commute with each other. The authors of [Tl] consider the matrix T{—u) with the entries 



-u 



-u + X 



hhijt 



Representing this matrix in the form 

T{-u) = e{-u)e "=0 



'Lo{u;X)e '==0 



(2.11) 



we obtain a dynamical elliptic L-operator Lq{u; X) over the algebra 1 = eT-,;i(0[„)[[c?A]] with 
Cartan elements h = (/ii, . . . , /i„), where C[[9a]] = C[[c?Ai , • • • j 9a„]], the elements dx,^ = 
commute with hi and do not commute with t,- 



2.2 Manin matrices and L-operators 

The RLL relations allow to construct the Manin matrices and g-Manin matrices investigated 
in [U m [15] starting from the L-operators. In particular, the dynamical RLL relation with 
Felder R matrix leads to the Manin matrices (here q = 1). We use the properties of these 
matrices to prove the commutativity of some function family tm{u) which will be constructed 
in the next subsection. We shall suppose in what follows that all matrix entries belong to some 
non-commutative ring which satisfies basically the conditions described precisely in [5l l6| [T5]. 

Definition 1. An nxn matrix M over some (in general, non-commutative) ring is called Manin 
matrix if the elements a, b, c, d of its any 2x2 submatrix, where 



M 



a 


. b ... 


c 


. d ... 



satisfy the relation ad — da = ch 
each other: ac = ca, bd = db. 



J 

be and their elements from the same column commute with 



Let Sm be the symmetric group and vr: 5„ — > End(C")'^™ be its standard representation: 
T:{a){vi (8) ■ ■ ■ (8) Vm) = ^0—1 (1) • • • <8> ^o— i(m)- Let us introduce the operator 



where (—1)'^ is a sign of the permutation a. We also use the notations A!^^'' 
(Am_fc)('^+^'---'™). For example, 



2^ 



[Eii (g) Ejj - Eij (g) Eji) . 



(2.12) 



In terms of the matrix (j2.12p Definition [T] can be reformulated as follows. The matrix M is 
a Manin matrix if and only if it satisfies the relation 



^(12) ^j(l) JV,/{2) = ^(12) m(^)m(2)^(12). 



(2.13) 
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Proposition 1. If L{u; X) is a dynamical L-operator then the matrix 

M = e'^^^L{u] A)e^^ = LD{u)e^^ , (2.14) 



where D\ = ^ E^k^^, is a Manin matrix. 
k=i 



Proof. Substituting u = —h to (j2.2p we obtain the formula 

R{-h; A) = V {Eu Ej, - E,j ® Eji) , (2.15) 

from which we see R{-h;X) = R{-h- X)A'^^^\ Consider the relation (f2lil]) at n — f = —h. The 
multiplication of this relation by A^^'^^ from the right does not change its right hand side, hence 
this does not change its left hand side: 

A)^a) (^)^g) ^u + h)= R^^^\-h; A)lS^^ (n)Lg) (n + /i) A^^^) _ (2.I6) 

Multiplying (j2.16p by e'^^a^ from the right we obtain 

A)mWm(2) = i?(i2)(_;i; A)M«M(2)^a2). (2.17) 

Let us note that the matrix (j2.15p is r elated to A^^^) by the formula B{X)R{-h;X) = yl^^^)^ 
where 



2^^(A, + n)^^^^- 

So, multiplying (pTT]) by B{X) from the left one obtains (pJ3]) . ■ 

Lemma 2 ( [5l [6] ) . If M is a Manin matrix invertible from the left and from the right then its 
inverse is also a Manin matrix. 

In particular, the matrix inverse to (j2.14p having the form 

{e-^^^L{u; X)e^^)~^ = e'^^Liu; X)~^e^^^ (2.18) 
is a Manin matrix. 

Lemma 3. If M is a Manin matrix then it satisfies the relations 

^[m,Af]^j(m+l) . . . j^.^{N) ^ ^[m,N] -^{m+l) . . . (A^) ^[m,Af] ^ ^2.19) 
j^[m,N\jy^{N) . . . ^ ^[m,JV]^(Af) . . . ]^{ra+l) ^[m,N] ^ ^2.20) 

where m < N. 

Proof. The idea of the proof of (12.19P is the following. It is sufficient to prove that the left 
hand side does not change if we multiply it by (— l)°''''+i7r((7j^j+i) from the right, where (Ji^i+i is 
an elementary permutation. But the last fact follows from (j2.13p . (See details in |15j.) 

Multiplying (j2.19p by (— l)°"'''*7r((T/st) from the both sides, where aist is a longest permutation, 
one yields the relation (|2.20p . ■ 
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Corollary 1. If L{u;X) is a dynamical elliptic L-operator then it satisfies 
^[m,N] Yl L'^^)(u + h{i-m-iy,X + hJ2EA 

m+l^i^N \ l=i+l / 

= aI'"'^] yI l(^) L + h{i-m-iy,X + hJ2 ^^""'""^ (2-21) 

m+l^i^N V l=i+l J 

N 

\ l=j+l 

= A^m,N] Yl L'^^Uv + h{j-m-l);X + hJ2E'^'A (2-22) 

N^j^m+l \ l=j+l J 

Proof. The relations (j2.2ip follows from the formula (j2.19p for Manin matrix (j2.14p : M = 
e-^^^L{u-X)e^-s^. The relations ^(TM follows from ([2:20]) for Manin matrix (121^ : M = 
e-^^L{u-X)-^e^^K ■ 

Let us consider the following matrix 

M['-'^]({7xa; {v,}; A) = n ^^'^ ( i^^}! ^ + ^ E ^^'1 =11 11 (2-23) 

l^i^m \ l=i+l / Is^i^m 

where 

(m N 
i=i+i i=j+i 

Taking into account ()2.5p we obtain the relations R^'^^ R'^"''^ = R^''^'iR^'^\ where 1 ^ i < a ^ m, 
m + l^j<b^N. Using it we derive the formula 

m[™'^]({^,};{z;,};A) = f[ R<^^^ I {u,};vf, X + h ^ 

N^j^m+l \ l=j+l 

= n n ^^-^^^ 

Let us set = u + h{i — 1), f j = f + — m — 1). Then, substituting M('^)(n; {wj}; A) and 
M(°)({^ii};^'; A)~^ as L-operators to the formulae ^^1711} (with m = 0, N = m) and ([2:^ 
respectively, taking into account the equalities (j2.23p and (|2.24p we obtain 

^[0,m,] j^[m„7V] ^. ^ ^[0,™] j^[m„7V] ^. ^)^[0,m] ^ ^2.25) 
j^[m,N]^[m,N] ^. ^) ^ ^[m,JV] j^[m,,7V] ^. ^)^[m,7V] ^ ^2.26) 

where rI'"'^] (n; t;; A) = = u + - = + - m - 1)};A). Similarly 

substituting M(o)({'Ui}; I); A)-^ and M(°)('u; {vj }; A) to the formulae (fOT]) and ([232]) (with m = 0, 
= m) respectively, we obtain 

^[o,m]M[™,7V](^.^.^)-i =^[0H]e^M(^.^;.a)-1aM^ (2.27) 

^[m,Ar]]^[m,7V](^.^.^^-i = ^[™,a^]]r[™,a^] A)-^ ^l"^'^] . (2.28) 
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2.3 Commutative families 

The integrals of motion for an integrable system related with an L-operator are obtained often as 
coefficients of a decomposition for some operator functions. These functions are constructed as 
traces of tensor products of L-operators multiplied by the "alternator" A''^'™] . In the dynamical 
case one should consider the tensor product of the operators (|2.9p . Now we use the facts proved 
below to establish the commutativity of these operator functions. 

Let us first fix the dynamical elliptic L-operator with Cartan elements /i^ and assign the 
notation L{u; A) to this L-operator. Its values at points u belong to the algebra EndC" (K> 

Introduce the following operators 

where m < N . Using the dynamical RLL relations (j2.7p and taking into account (j2.6p we derive 
the commutation relations for them: 

/ ^ \ 



{u,]-{v,]-\-nY,E^'^ LM({na;A)L['"'^l({^;,};A) 



\ 1=1 / 

= A)LM({ni}; A)M[™'^] ({u^}; {vj}; A + hh). (2.29) 

Setting Ui = u + h{i — 1), vj = v + h{j — m — 1) in this relation we obtain the main relation 
which we need to construct a commutative family: 

M[™'^1 (^u;v;X-hJ2E^^^^ LM(n; A)L[™'^1(i;; A) 
= lI'"'^! {v; A)l[°'™] (n; A)r[™'^1 {u;v;X + hh) , 

where we set L['''^](n; A) = L['''^]({ni = u + h{i - a - 1)}; A) for a < b. 

Let A„ = C((Ai, . . . , A„,)) be a completed space of functions. The operators Dx acts on 
the space A„ (g) C"-, so that the operators L["'''1(u; A) act from the space A„ (g) (C")®(^~") to 
A„ (g) (C'')(S>(b-a) ^9^. if It is a fixed point then h^^'^^u; X) G End(C")®(''-") 21^, where 
2ln = A„[e^'^^] (g) yi. Consider the subalgebra f) C 9^ C 21^ generated by elements hk and its 
normalizer in 21^: 

Let us note that 2l„[) is a two-side ideal in 
Theorem 1. Let us define the ^n-valued functions 
i^(t/) = tr(AML[0H(^.A)), 

where the trace is implied over all m spaces C". They commute with Cartan elements h^: 

hktm{u) = tm{u)hk. (2.30) 

Hence they take values in the subalgebra And they pair-wise commute modulo the ideal 

2l„f) C mn.: 

tm{u)ts{v) = ts{v)tm{u) mod 2lnf). (2.31) 
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Proof. The relation ()2.30p follows from the formulae 

m 

[/i„LM(n;A)] = - J][4),LM(n;A)], 
1=1 

m m 

S« AP'H = ^[0.™] J2 Efj (2.32) 
1=1 1=1 

and from the trace periodicity. 

Let us stress that the commutativity (|2.31|) modulo 2l„f) is a corollary of the formula 

tr (a[0.-U[-,a^]l[0H(^. A)L["^'^l(t;; A)) 

= tr (^[0,H^['n,Af]jL[m,7V] x)h[0M a)m[™'^1 (n; v; X + hh)R^'^'^ (n; v; A)"^) , (2.33) 

where N = m + s and the trace is considered over all N spaces C". This formula can be deduced 
as follows. Consider the left hand side of ()2.33p . Applying the relation (|2.29p we can rewrite it 
in the form 

/ / iv \ -1 



X lIO'*"] (n; X)R^"''^^ (n; v; X + hh)j . (2.34) 

Note that by virtue of the formula (|2.32p the relations (j2.27p , (j2.28p with the shifted dynamical 

N 

parameters A — > X—h ^ are still valid. Therefore, we can insert the idempotent A^'''™]^!"*'^! 
1=1 

, N s-l 

after R[™'^1 (u;v;X-hY, to the expression I^Ml- 



1=1 ' 

/ TV \ -1 



1=1 



X lP'^I (n; A)M["*'^1 (n; v;X + hh)j. 

Similarly, taking into account (|2.2ip we can insert A[°'™U["*'^1 after L[°''"1(u; A) and, then, usmg 
a^TM . a^TM and the trace periodicity we can eliminate ^["''"UI™'^] before R^""^^^ (u;v; X - 

n • Thus we obtain 



1=1 



tr rI'"'^] I u 



N 

■v-X-hY.E'^'^ ^[0'-UMLK^](t;;A) 
1=1 / 

X LM(n; A)A[°''"U['"'^]m[™'^](u;i;; A + hh)^ , (2.35) 
where one can delete after L['''"^l(n; A). By virtue of (|2.32p we can present the part 



JV 



after rI™-'^] (u;v;X — hY E^^'^ J as e '=i ^('^)> where the matrix Y{X) does not contain 



=1 
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the shift operators. Using the formula 

tr X(A)e '=1 'y(A) =tr e '=i ^ YiX)X ix + h^E^^'^ 

where the matrix X{X) does not contain the shift operators and elements of the ring IH, one can 
rewrite (|2.35p as the right hand side of (j2.33p . ■ 

Thus we obtain the family of functions tm{u) with values in the algebra Their images 
by the canonical homomorphism 97„ — > DT^/Slrtf) commute and we obtain the commutative 
family of functions in the algebra OT^/Slnfl- Decomposing these functions in some basis functions 
(pmjiu) one yields a commutative family of elements Im,j of this algebra: tm{u) = Im,j4'k,jiu) 
mod 2l„j). In this way one can construct an integrable system corresponding to a given L- 
operator. 




2.4 Characteristic polynomial 

The functions tm{u) can be gathered to a generating function called the quantum characteristic 
polynomial. This polynomial is defined in terms of determinants of Manin matrices. Considering 
the Gaudin degeneration we shall calculate the degeneration of the characteristic polynomial, 
although the degeneration of the functions tm{u) can not be obtained explicitly. 

Introduce the following notion of a determinant for a matrix with non-commutative entries. 
Let M be an arbitrary n x n matrix, define its determinant by the formula 

det M = (-l)"M,(i),iM,(2),2 • • • M,(„),,. 

(r£S„ 

Usually this determinant is called the column determinant. One can show |15j that for a Manin 
matrix M its column determinant coincides with the completely symmetrized determinant and 
can be represented in the form 

det M = tr (a[°'"1m(^)m(2) . . . M^) . (2.36) 

Proposition 2. Let us consider the Manin matrix M = e~^'^^L(u; A)e^'9^. The matrix 1 — M 
is also a Manin matrix and its determinant generates tm{u) o,s follows: 

^ n 

P{u,e^^) = Aet {I - e-^'^^L{u;X)e^-s^) = ^ (-irt„('u)e"''^li , (2.37) 

m=0 

where tQ{u) = 1. This implies that the determinant (j2.37p commutes with hk as well as that it 
commutes modulo 2l„f) with itself for different values of the spectral parameter: 

[P(u,e'^"),/ifc] = 0, [P{u,e^^),P{v,e^-)]=Q mod2t„[). 

Proof. First we note that the functions tm{u) are related to the principal minors of the mat- 
rix M as follows 

det m;;-;;^ = tri,...,^ (^a^^Mm^^) ■ ■ ■ M^"^)) = tmiu)e'^^^, (2.38) 

l^il<---<im^n 

where M^^'''^'^ is the corresponding submatrix. Due to the formula (|2.36p the left hand side 
of ()2.37p decomposes as follows 

n 

det(l -M) = Y i-^r Yl *^i.-..,n (^["'"iM^'^i^M^^^) . . . m(^™)) . (2.39) 

m=0 l^ki<-<kmiin 
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Using the recursive relation = - (n - l)^[°''"-i]7r(cr„_i,„)A[0'"-il) one yields 

tri,...,n (A[o."1m('^i)m('^2) . . . jvf (fc™)) = (A[Mm(^)m(2) . . . M^™)) (2.40) 

n! ' ' 

The formulae (|238]) . (|239]) . ([TIO]) imply (|237D . ■ 

2.5 Newton identities and quantum powers 

In the theory of classical integrable systems commutative families are usually provided by traces 
of powers of a classical L-operator. Their quantization can be presented as traces of deformed 
powers of the corresponding quantum L-operators. These deformed powers are called quantum 
powers of L-operator. For rational L-operators these quantum powers were described in details 
in [5]. The main tool to obtain the quantum powers are Newton identities for Manin matri- 
ces. They allow to express the quantum powers of L-operators through the coefficients of the 
characteristic polynomial. 

Lemma 4 ( [3 E] ) . Let M he a Manin matrix. Introduce the elements 

Qm = tri„„,™ (a[o,™.]^(i) . . . = det M;i;-^, (2.41) 

go = 1 and qm = for m > n. Then for any m ^ the following identities hold 

m—l 

mqm = Y^i-ir^'^'qktv {M^-^). 

k=0 

They called Newton identities for the Manin matrix M . 

Consider the Manin matrix M = Lj:,{u)e^'^ = e~^^^L{u]X)e^'^ . The elements (I2.41j) for 
this matrix have the form q^n = tm{u)e"^^~a^: see the formula (12.380 . The powers of M are 

= L'§{u)e^^^, where h'^iu) = Ld{u)Ld{u + h) ■ ■ ■ Ld{u + {k - l)h). 

The matrices L^'(n) are called quantum powers of the L-operator L/)(n). The Newton identities 
allow recursively express the traces of the quantum powers tr (L^^(ti)) via the functions tm{u) 
and vice versa. Thus we obtain another commutative family of functions generating the same 
commutative subalgebra in 9T„/2t„f): 

[tr(LM(n)),tr(Lg^l(7;))] =0 mod 2l„f). 

2.6 Trigonometric limit 

Up to now we have considered -L-operators satisfying dynamical i?LL-relations with a certain 
dynamical elliptic i?-matrix. But the presented technique is sufficiently universal to be applied 
to many other cases. This approach is directly generalized to the case of arbitrary dynamical 
i^-matrix such that R{-h; A) = ^(A)^^^^) fo^. 

some invertible matrix B{X). The corresponding 
L-operators give us a family tm{u) commuting modulo %rS)- In particular, it works for the 
trigonometric and rational degenerations of the Felder i?-matrix. Using the same scheme for 
their non-dynamical limits we obtain commuting functions tmiu) (modulo 0). 
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Here we briefly discuss the dynamical and non-dynamical trigonometric limits. The Felder 
i?-matrix in the limit r — > ioo takes the form 

-1 " 

^ , , , zq — wq \ - ^ ^ 

-Rtrigl'", V, A) = } Eii(g) Eii 

z — w ^-^ 

i=l 

+ E ( ^^^^g.. « E„ + '7 "'';<;- E,.) . (2.42) 

where z = e^'^*", -u; = e^'^*'', q = e'"^^, fikj = e^'^'^^'^K An L-operator satisfying the relations (12. 7p 
and (|2.8p with this i?-matrix - a dynamical trigonometric L-operator - defines the commutative 
family tm{u) (modulo corresponding 2tnf)) by means of the formula (j2.37p . for instance. One can 
also consider the commuting traces of quantum powers in the same way as in the Subsection 12. 5[ 
Consider the limit Xk — Xk+i — > —ioo, that is fiij oo for i < j and fiij — > for i > j. The 
limit of the matrix (j2.42p is a non-dynamical trigonometric i?-matrix 

_ _i n / 
^trig('W -v) = Ruigiz, w) = ^ Eii ® -Eji + ^ I qEii (g) Ejj + q^^Ejj (g) En 

1=1 i<j \ 

+ ^-^f^E, ^ E,. + E.,] . (2.43) 

[z — w) [z — w) I 

Let L{z) be an L-operator satisfying the usual i2LL-relations with this matrix: 

Rtrigiz,w)L^^\z)L^^\w) = L^^Hw)L^^Hz)Rtrisiz,w). 



The matrix M = L{z)q is a Manin matrix. The corresponding characteristic polynomial 
det(l - L{u)q^^~^) = Yl (-l)*m(-2)9^'"^s^ defines a commutative family: [tmiz),tm{w)] = 0. 

■m=0 

Then, the quantum powers in this case have the form LM(z) = L{z)L{zq^) ■ ■ ■ L{zq^^^^'^'>), 
[tr (LM(2)),tr (lH(u,))] =0. 
Let us consider the tensor 

i?trig(^, w) = F(2i)i?t,ig(z, w)F~^ = ^'^ L;ii ®Eii + ^[ Eii L^ij 

2 = 1 2<J \ 



{q^^q_^ {q-q ^)z 

^Eij ® Eji H - 

[z — w) [z — w) 



+ Ejj ® Eii + ^\ Eij Eji + ; ® Sij- , (2.44) 



where 



F = q =2_^Eii(^ Eii + 2^(q^Eii(^ Ejj + q 2 Ejj O Ei 

i=l i<j 

This is a standard i?-matrix for the quasi-triangular Hopf algebra f/g (§[„). This algebra can 
be described by two L-operators L^{z) G EndC" (g Uq{Q{^) satisfying i?LL-relation with the 
i2-matrix (j2.44p . Each L-operator L{z) G EndC" 9^ satisfying i?LL-relation with this R- 
matrix defines a homomorphism from the certain Hopf subalgebra ^ C f7q(gl„) to 9^H The 
i?-matrix (j2.44p considered as an L-operator defines a representation vr^ : *B — > EndC^. The 



The algebra 23 is the algebra described by one of those L-operators, for example, by L'^{z). In this case the 



homomorphism is defined by the formula LfJz) 1— » Lij{z) 
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subalgebra *B contains elements hi, . . . ,hn such that '/r^(/ifc) = Ekk, ^{hk) = 0, A(/ifc) = hk ^ 
1 + 1 ® /ifc, where e and A are co- unity and co- multiplication of The element 

T = q 

satisfies the cocycle condition (the Drinfeld equation) (A (g) id)(J^) = J^'^^'^^ (id (8)A)(:F), the 
condition {e®\d)J- = (id ®£)J- = 1 and, therefore, defines a twist between two co-multiplications 
of the algebra C/g(0[„), the standard one A and twisted one A^ defined as l^^{x) = J^A(x)J^-^ 
Let (f)-. 55 ^ $H is the homomorphism defined by a given L-operator L{z). Define the matrix 

G = q =l^q ^^='+' ^ Eii. 

1=1 

where /i^ = (f>{hi^). From the theory of quasi-triangular Hopf algebras it follows that the twisted 
L-operator L{z) = GL[z)G satisfies i?LL-relations with the matrix (j2.43p . This means that each 
L-operator L[z) defines the commutative family tm.{u) by means of the characteristic polynomial 
for the Manin matrix M = GL{z)Gq^^s^ . 

Let us also remark that the commutative families for L[z) can be constructed without twist- 
ing. In this case one should use some generalization of the notion of Manin matrices. The 
matrix M = L{z)(f'^'s^ belongs to such kind of generalization called q-Manin matrices (see [15] 
for details). Nevertheless the approach of g-Manin matrices is a little more complicated. Let 
us also remark that the same situation occurs for the C^5,p(0l„) L-operators considered in [16], 
there exist a dynamical elliptic analogue of the twist J- relating the C/g^p(0l„) i?-matrix with 
the Felder i?-matrix. So we expect that there exists a generalization of the g-Manin matrices 
corresponding to the dynamical elliptic quantum group C/g,p(0t„)Il 



3 Characteristic polynomial for elliptic Gaudin model 

We consider a degeneration of the dynamical elliptic RLL relations: ?i — > 0. In particular, this 
degeneration describes the dynamical elliptic 0l,„ Gaudin model. To relate this degenerated L- 
operator with the elliptic Gaudin L-operator in the standard formulation we introduce a twist of 
this L-operator. Degenerating the commutative family obtained in the previous section we yield 
a commutative family for the elliptic Gaudin model. We see that the obtained results generalize 
the elliptic 5(2 Gaudin model case investigated in [TSl [T9] . 



3.1 h — > degeneration 

Here we consider the degeneration of the previous section formulae obtained in the limit — > 0: 
we compare the coefficients at the minimal degree of h that does not lead to a trivial identity. 
In the degeneration of the dynamical i?LL-relation, for example, we consider h?. And the 
determinant (j2.37p is proportional to the n-th power of h. We avoid using the term the classical 
degeneration because we don't want to mismatch the considering degeneration with the classical 
mechanics degeneration. 

Suppose that L(n; A) is a dynamical elliptic L-operator of the form 

L{u;X) = l + hC{u;X)+o{h). (3.1) 

^In some formulation Uq,p{gl„) should be interpreted as a Hopf algebroid [17] . 
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where the entries of C{u; A) belongs to the algebra 5Ho = The matrix C{u; A) is called 

dynamical classical elliptic L-operator. It satisfies the dynamical classical rLL-relations 

[£(1) (n; A) - 5« , £(2) x) - 5f ] - /i,|-r(n - ^; A) 

fc=i 

= [£W(n;A) +/:(2)(?;;A),r(u-v;A)] (3.2) 
with the dynamical classical elliptic r-matrix 

^(-; ^) = ^ E ^ + E (^^- ^ + ,(;),(-4) ^^ ^ J • (^-^^ 

i — 1 iy^j 

The matrix (13. 3|) is related with the Felder i2-matrix ()2.2p via the formula 

R{u; A) = 1 + hr{u; Xn) + o{H). 
The relation ()3.2p follows directly from the relation (j2.7p . In turn, the relation p.Sp implies 

(Eii + /ii)/:(n; A) = /:(n; A)(Eii + hi). (3.4) 



Let us remark that two different Lie bialgebras described by the relations (13. 2p and ()3.4p are 
considered in great details for the s[2-case in [20]. 



Lemma 5. If a Manin matrix has a form M = 1 + hAi + o{h) then Ai is also a Manin matrix. 

Proof. Substituting the decomposition of M to (12.13P we obtain the equation (I2.13P for M. at 
second power of h. ■ 

Applying this lemma to the Manin matrix M = e~^^^L{u; X)e^d^ we conclude that the 
matrix -M = — D\ + C.{u; A) is a Manin matrix. Theorem [1] and Proposition [2] give us the 
following theorem. 

Theorem 2. Let An = '^q ® A.„[(?a] and Nn = ^^„(f)) = {x G An \ C An^}, where 
An = C((Ai, . . . , An)). Define the Mn-valued functions Sm{u) by the formula 

— -Dx + C{u;X)j = Y,Sm{u)i — ] , (3.5) 

^ m=0 ^ 

where so{u) = 1. They commute with Cartan elements hj-: 

hkSm{u) = Sm{u)hk (3.6) 
and they pair-wise commute modulo An^: 

Smiu)si{v) = si{v)smiu) mod Anh- (3.7) 

The decomposition coefficients of the functions si{u), S2{u), . . . ,Sn{u) form a commutative 
family in the algebra 7V„ at the level /i^ = 0. This means that their images by the canonical 
homomorphism 7V„ Nn/ An^ commute with each other. Below we consider some special class 
of L-operators C{u; A). 
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3.2 Elliptic Gaudin model 

Theorem [2] provides the functions Sm{u) generating the quantum integrals of motion for the el- 
liptic Gaudin model. This model is defined by a representation of the loop Lie algebra z~^]. 
Let us first consider a homomorphism p: U{Qln[z, z~^]) — 9\q from U{gln[z, z~^]) to some al- 
gebra 5Ho- To this homomorphism we associate the elliptic half-currents 

where i ^ j, {cj-,} is a basis in gl^ for which Cij i— > Eij is a faithful representation. The set 
{ciiZ^} is a basis in Qln[z, z~^]. The matrix C{u; A) with the elements 

C,,iu;\)=e±iu;X), A) = (t,; A) + V ^^/i^, (3.10) 

where = p{ekk), satisfies the relations ()3.2I) and ()3.4p . 

Consider an important example of the homomorphism p. The evaluation homomorphism 
Pv'- C/(0[„[z, z~^]) U{q{j^) defined as 

Pv{(^ijZ ) — V eij, 

where i,j = 1, . . . ,n and f G C is a fixed complex number. Using the standard co-multiplication 
of the universal enveloping algebra [/(gl^Jz, z~^]) one can construct the more general evaluation 
homomorphism p^^j : U {glni^i ^~^]) ~^ ^ (gt^)*^^ in N fixed points vi, . . . , u at G C: 



k=l 



N 



where e\j^ = ^ eij ® i'S{N-k) -g ^^iq basis element corresponding to the A;-th tensor 

component of C/(0[„). For the given representations 11^: [/(fll^) — > EndV^ one can construct 
the representation H^^}- U{Qln[z, z~^]) —>■ EndV, where n|^} = (IIi ® ••• (g) Hat) o pj^j. and 
V = Vi iSi • • • (d) Vn ■ In the representation ^{v} the half currents p.Sp , (|3.9p take the forms 

N an _ \ 

4(-;A) =e^(n) =5:^)j^n,(e..)W, (3.11) 



(9(n - Vk) 
9{u -Vk-\- Xi 



k=l 

N 

4(M)^Ei:^n.(e..)'". (3-) 



where Ilk{eij)'^^^ = l^C^-^) Ilk{eij) ® i<^{N-k) ^^j^g tensor component). The matrix 
C{u; A) defined by (j3.10p with the currents ()3.1ip . (j3.12p is called the elliptic Gaudin L-operator 
corresponding to the representation n|^}. 

In the case of homomorphism p = pj^} the corresponding commutative family can be ob- 
tained as follows. Let Lq{u;\) be the dynamical elliptic L-operator defined by (12. lip . Shifting 
the argument we again get a dynamical elliptic L-operator Lq{u — v;X). This is also a ma- 
trix over T = eT-,?i(0(n)[[c^A]]- Let A„ = A„[[^]] [[5a]] C T. Let Lk{u — Vk] A) be the matrices over 

= T . . . (^^^ 1 with the entries {Lk{u-Vk; X))^. = l®('=-^)o(Lo(u-t;fc; A)) .^.(g)l®(^-^) and 



N 
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h'' = (h'l,. . . ,h^) withh'^ = 1'^^'' (g)/ij(g)l®(^ Then applying Lemma[I]we obtain a dynam- 

ical elliptic L-operator (n; A) = Y[ Lk{u — Vk]\ + h ^ /iM over d\ with the Cartan ele- 

N^k^l ^ l=k+l ' 

N 

ments /i = Representing L^j\f^{u; A) in the form p.ip we obtain an L-operator jC[^](n; A). 

k=l _ 

This is a matrix over 9\q = To (8)^ • • • ®^ To, where Tq = %/h% and = A,„[[(?a]] C Tq. Due 

N 

to Theorem [2] the functions s"m\u) defined by (|3.5p with this L-operator C^]^t^{u;\) satisfy (|3.6p 
and (j3.7p . Let v?: 9lo — > So be a homomorphism to an algebra So- Applying it to the equa- 
tions (j3.5p . ()3.6p and ()3.7p we conclude that the functions Sm{u) = ip{sm\u)) coincide with 
the functions defined by (j3.5p with the L-operator >C<^(n;A) = (/^(^[tv] ("w; A)) and satisfy (|3.6p 
and (|377ll . 

Consider the homomorphism tpQ-. 1q ^ (/{qI^) A.„ defined by the formulae 
'^ofe) = E ^^'^'^ - ^' = ^e., for z ^ j, 

^o{hk) = Cfcfc and (po{dx^,) = d\^. Let : IHo ^ ^^(s^n)®^ be a homomorphism defined 

as 93[Ar] = V'o <X';4 ■ ■ ■ '^A ''^o- Then the operator C^^j^^{u;\) = Lp\^N]{C[N]{'^\ ^)) coincides with 
^ — ' 

N 

the L-operator (j3.10p corresponding to the homomorphism p = The functions Sm""^ {u) 

defined by (j3.5p with £1^^^^ (u; A) satisfy (j3.6p and (j3.7p . Analogously, considering 93 = (Hi (g) 
• • • (8) Htv) o (^[at] we obtain the Gaudin L-operator corresponding to the representation = 
(Hi (8" • • • Hat) o /3|^}. Thus we derive the following proposition. 

Proposition 3. Let Ilk - U{Qln) ^ End 14 be representations andC{u;X) be the elliptic Gaudin 
L-operator corresponding to the homomorphism Hj^j = (Hi (8 • • • (8 Htv) o P{v}- Then Sm{u) 
defined by the formula (I3.5P with this C{u; X) satisfy (13. 6p and (I3.7p . These are integrals of 
motion for the dynamical elliptic qI^ Gaudin model. 

So we have shown that the L-operator (j3.10p corresponding to the certain representation 
p: U {Qln[z , z~^]) —>■ fHo = EndF defines a commutative family Sm{u). We conjecture that this 
is true for any homomorphism p. 

Conjecture 1. For any homomorphism p the corresponding L-operator (j3.10p defines a com- 
mutative family Sm{u) by the formula ()3.5p . that is Sm{u) satisfy (I3.6p and (I3.7p . 

For the case n = 2 this conjecture can be proved using a certain dynamical elliptic L-operator 
constructed in [21]. In notations of that paper this is the L-operator L(tt,A) constructed from 
the L-operator L'^{u). The entries of L{u,X) belong to the certain algebra 9^ = C/a5 (via some 
embedding) and the degeneration of this L-operator gives the L-operator (j3.10p for p = id|l| 
Thus the L-operator for p = id defines functions s}^{u) satisfying ()3.6p and (13. 7p . The L- 
operator for the arbitrary homomorphism p defines function Sm{u) = p[sm.{u)), applying p to 
the relations (j3.6p and (j3.7p for s5^(n) we conclude that Sm{u) satisfy the same relations. 

For the general n this conjecture can be proved in the same way. For example, it follows from 
the Conjecture 5.1 of the paper [16]. There the authors present an analogue of the L-operator 
L^{u) for the quantum group f/g,p(s[„). The existence of a twist relating this L-operator with 
an analogue of L{u, A) can be easily proved. 



*The degeneration of this algebra $H = UnQ with zero central charge K = and without co-central charge is 
SHo = $H/ftlR = (7(£|[„[2,2-i]). 
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Theorem [2] with Proposition [3] provide us a method to construct the quantum elliptic Gaudin 
model for general n in analogy to the rational case treated in [TJ]. Let us firstly note that the 
connection whose determinant Q{u, du) provides a generating function for quantum Hamilto- 
nians is the KZB connection. Moreover this construction allows to establish important relations 
of the elliptic Gaudin model with the Langlands correspondence program. The construction of 
the characteristic polynomial Q{u,du) called the Universal G-oper in this case provides us 

• a universal construction of the commutative family; 

• a universal description of the eigenvalue problem; 

• a relation between the solutions of the KZB equation and the wave-functions of the Gaudin 
model; 

• a relation with the centre of the affine algebra at the critical level 

analogously to the results of [22\ [23] in the rational case. 

Usually the elliptic Gaudin model defined by the L-operator 

n 
ij=l 

with the currents (jS.lip . (j3.12|) or, more generally, with the currents (|3.8p . (j3.9p . In the next 
subsection we describe the relation of the L-operator (I3.13P with the L-operator C{u; A). 



3.3 Twisting of a classical dynamical L-operator 

First we consider the twists of dynamical r-matrices for arbitrary Lie algebras. Let g be a (in 
general, infinite-dimensional) Lie algebra, g (8) be a tensor product completed in some way and 
f) C g be an n-dimensional commutative subalgebra called Cartan subalgebra. Let us fix a basis 
{hi, . . . , hn} of [). An element r(A) G g (8" g depending on TL dynamical parameters Ai, ■ ■ ■ , A^ is 
called a dynamical classical r -matrix if r^^^^ (A) + r^^^^ (A) G (g (X" g) ^ and it satisfies the classical 
dynamical Yang-Baxter equation 

[[r(A),t(A)]]+P,(t(A))=0 (3.14) 

where 

[[a(A),6(A)]] = [a(i2)(A),6(i3)(A)] + [a^'') (X) , bi''\X)] + [a('^\X),b('^HX)], (3.15) 



The space (g (8> g) ^ consists of the elements x G g g such that [x, y 1 + 1 y] = for all y G g. 
The brackets in the right hand side of (|3.15p mean the commutators in U{q) ® C/(g) (8) f^(g)- 

Lemma 6. Let r(A) and f(A) be elements o/g g. Suppose that f'^^^^(A) + f^'^^\X) G (g g)^ 
and 

[[f(A),f(A)]] + [[?(A),f(A)]] + mX),i{X)]]+V,{f{X)) = 0. 

Then, r(A) is a dynamical classical r-matrix if and only if the matrix r(A) = r(A) + f(A) is 
a dynamical classical r-matrix. The element f(A) is called classical dynamical twist. 
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Consider the loop algebra g = Qln[z, z^^\ with Cartan subalgebra f) spanned by the elements 
^fc = ^kk- This is the Lie algebra of constant diagonal matrices. Let r(A) be a dynamical 
classical r-matrix and let ■Ku'- ^^n\Z-,z~^\ EndC" be the standard evaluation representation 
defined as 7r„ : eijZ^ i— EijU^ . Since vr^ is faithful the matrix r can be represented by the matrix 
r(n, v; A) = (vr^ (g) ■Ky)x:{\) G EndC" ® End , v, v~^]], where the dependency of u and v 
is understood in the formal way. Each homomorphism p: U {gl^lz , z~^]) — > 9^o defines the L- 
operator C{u;X) = (vr^ (S) /o)t(A) satisfying the dynamical rLL-relations ()3.2p with the r-matrix 
r{u,v;X). 

Let r{u,v; A) S EndC" ® EndC"[[f]]((ii)) be the r-matrix r{u — v; A) (formula (jS.Sp ) under- 
stood as a formal series with coefficients in End C"® End C". Relation (j2.3p implies that the cor- 
responding element r(A) is a dynamical classical r-matrix. The matrix C{u; A) = {tTu^ p)t{X) co- 
incides with the L-operator (j3.10p corresponding to the homomorphism p: U {sini^, z~^]) 9^o- 

Consider the following twist 

It satisfies the conditions of the Lemma [6l Indeed, the twisted r-matrix r(A) = r(A) — f (A) is 
represented by the matrix 

0^(11) ^ 0(11 A''l 

fiu, V- A) = r(n, vX) - (vr^ ® 7r^)f (A) = V En <^Eu + y] f Eij ® Eji, 

ti{u) ^ ^ 0[u)0[-Xij) 

where (7ru07r^)f(A) = ^ -^rj^Eii® Ejj is the representation of the twist (j3.16p . One can check 

the equalities [[f(A), f(A)]] = [[r(A), f(A)]] + [[f(A), r(A)]] = Vx{^{X)) = applying homomorphism 
(tTu ® 7r„). So that the matrix r(A) also satisfies (|3.14p . 

The twisted L-operator C{u; A) = (vr^j ® /o)t(A) satisfies the rLL-relations (13. 2p with r- 
matrix f{u,v]X) and the relation (j3.4p with = p{hk) = p{ekk)- It coincides with the L- 
operator (13.13P related with the L-operator C{u; A) via the formulae 

Cij{u; A) = Cij{u; A), Cu^u; A) = £n(u; A) + -j^jr^hk, 
where j 7^ i. The characteristic polynomial from Theorem [2] in these terms takes the form 

Q{u, du) = det I |- - 5a + C{u; A) + E ^''J(j!^^^ I " ^^'^^^ 

Remark 2. Calculating determinant (I3.17P modulo ^nf)[[^i ""^j ^*]] in the n = 1 and n = 2 
cases we can omit the sum V En^aTr^h-i in the determinant. However for n ^ 3 one can not 

simplify it in this way, because even in the case n = 3 omitting this sum we lose the following 
term of (5(ti,9„) (mod u]] [9^]): 

'(A12) , (^'{^2z)\ p / x\r ( w 
+ Q,. N >Ci3(n; A)£3i(n; A). 



^(Al2) ^(A23) 

3.4 The SI2 elliptic Gaudin model 

Let us relate the generating function Q{u,du) in n = 2 case with the generating function 
describing the 3(2 elliptic Gaudin model [18l [191 [23] . The latter one corresponds to the QI2 elliptic 
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Gaudin model defined by a representation p: U{q12[z^ ^~^]) ~^ such that p: e^^^ (n) + (^) ^ 
0. For example, p = (p®- • •(8)p)op|,y|, where p(eii+e22) = 0, p(eii -622) = 611-622,^(612) = ei2, 
p{^2i) = 621- Then the L-operator takes the form 

Ziu-X)-( ^^^''^''^ 

where A = A12 = Ai — A2 and the half-currents are 

( \ +f \ +( \ \^0'{u-Vs) I (s) {s)\ 
h+{u) = e+ (n) - e+^iu) = gr^^ _ (^n " ^22)^ 

s=l 

e+(u) - e+ (u- A) - V + 

hN = e2ii-;>^) = Y. g(^u-vM->^f ^'- 

Since the L-operator depends only on the difference A = Ai — A2 we can restrict Q{u, du) to the 
space A = {a € A2 I ((9ai + dx^)a = 0} C A2 coinciding with C((Ai2)). Let A = "iJlo ^ A[(9a] 
then the values of the functions Sm{u) belong to J\f = ^^(f)) = {x ^ A \ hx ^ Ah}. Since 
p: /ii + /i2 — > the operator Dx in this representation has the form Hd\., where H = En — £'22- 
Now we calculate the column determinant (j3.17p : 

Q{u, du) = det (I; - 5a + An; A) - 

d.Q^ + h+{u)/2-'-^h/2 f^{u) 




m + dx-h+{u)/2-'-^h/2^ 

du) m d^~ '^''^' 

where h = hi — h2 and S\{u) is the following A/"- valued function 

Sx{u) = {dx - h+{u)/2f + duh+{u)/2 + e+{u)f+{u) mod Ah, 



which commutes with itself: [S'a(u), 5A(f )] = mod Ah. Using the commutation relation 

9(A) 



[e^(ii), /^(w)] = —-§:^h'^{u) + {^QT^)'h one obtains the generating function for the 5I2 elliptic 



Gaudin Hamiltonians: 

Sx{u) = {dx - h+{u)/2f + {et{u)f+{u) + /+(n)e+(n)) /2 mod Ah. 

This generating function was considered and used in [181 US] to find the eigenfunctions of this 
model. 

3.5 Quantum powers for the Gaudin model 

The quantum powers for the elliptic Gaudin models are defined in terms of the operator Cd{u) = 
C{u; A) — Dx. The first two quantum powers of this operator coincide with its ordinary powers: 
C^^\u) = 1, C^^\u) = Cd{u)- The others are defined by the following recursive formula 

^[k+i]. . ^ , s^[k]. X , dd^\u) 
Ojj \u)=Cd{u)OjJ{u) + — — — . 
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Considering the Newton identities for the Manin matrix = ^ + ^d{u) = ^ + ^{u; A) — Dx 
we conclude that the traces of the quantum powers form an alternative commuting family of 
functions for the elliptic Gaudin model: 

[tr(/:gl(^x)),tr(/:g^l(t-))] =0 mod Ani). 

The proof of this formula is the same as for the rational Gaudin model. See details in [5]. 

Instead of the quantum powers of Co{u) one can consider "simplified" quantum powers 
Cd{u)c}]^ '^\u). Their traces also form the commuting family of functions defining the same 
integrable system: 

[tT{CD{u)C^~^\u)),ti{CD{v)C^o~^\v))]={) mod An^- 

In particular, this family includes the trace of the ordinary second power 

i,{CD{uf)=t^{{C{u-\)-Dxf). 

4 Conclusion 

In Section [2] we have constructed a commutative family starting with the Felder iZ-matrix. All 
results of that section are valid if we would start with another dynamical iZ-matrix satisfying 
R{—h;X) = B{X)A^^'^^ for some invertible matrix B{X). The universality of this method also 
implies that one can consider an arbitrary L-operator satisfying dynamical i?LL-relations. In 
other words, it can be applied to a wide class of models which are described by using a dynamical 
i2-matrix. Moreover, the integrals of motion for the models corresponding to a non-dynamical 
iZ-matrix can be obtained in the same way, or more precisely as it was done in [13] for the 
rational case. In particular, it should work for the XyZ-model or other models described by 
the Belavin i?-matrix. 

The use of the Manin matrices is important to provide the universality of our method due to 
the universality of that notion. First, the characteristic polynomials of Manin matrices are key 
objects to obtain the commutative family of the (elliptic) gt^ Gaudin model by the degeneration. 
Then the Manin matrices are used to relate the constructed commutative families with another 
important class of commutative families - the traces of the quantum powers. The Manin matrices 
and their characteristic polynomials find their applications to many important problems in the 
theory of integrable systems. It is also convenient to use the properties of Manin matrices to 
prove commutativity of the constructed families. 

Analysing the trigonometric limit we have noticed that L-operators related with the Hopf 
algebra f7g(g[„), such as the L-operator of the XX Z-model, gives also a Manin matrix after 
some simple twisting. This means that the corresponding quantum determinants, characteristic 
polynomials, quantum powers and Newton identities can be written in terms of the ordinary 
anti-symmetrizers ^(^'•••''") by using the theory of Manin matrices, while consideration of the 
non-twisted L-operators leads to the using of g-deformations of these anti-symmetrizers and of 
a (^-analogue of the Manin matrices. 

In Section [3] we showed that the degeneration of each dynamical elliptic operator L(n; A) 
gives a commutative family for the corresponding degenerated L-operator C{u; A). To justify the 
obtained results for a given C{u; A) one have to prove the existence of the corresponding L(n; A). 
We considered the L-operators C{u; A) defined by some homomorphism p: C/(g[„[z, z~^]) — > 
or a representation p: U {gl^lz , z~^]) EndV. We discuss the existence of the corresponding 
L-operators L{u; A) presenting them explicitly or referring to other works. 

For the representations p = ^{v} we present the corresponding L{u;X) using the notion of 
the small elliptic quantum group suggested in [H]. To obtain the commutative family in the 



Manin Matrices, Quantum Elliptic Commutative Families and Elliptic Gaudin Model 



21 



case of arbitrary homomorphism p: C/(g[„[z, z^^]) it is sufficient to construct L{u; X) 

corresponding to p = id. In the (technically 5(2) case the L-operator L{u; A) corresponding 
to p = id in fact constructed in [21]. For the general gl„ (technically sin) case the (non- 
twisted) L-operator L{u] A) corresponding to p = id is written in [16] as a conjecture. One 
also need to relate the i?-matrix using in [16] with the Felder i?-matrix by some twist as in 
Subsection 12. 31 The existence this twist can be quite easily established using the ansatz = 
exp ^1 Si^^j ^ijWh'i (8) hj^ , where (pij{X) = —(pji{X) are C[[/i]] -valued functions. Let us finally 
remark that the degeneration of this twist gives the classical twist considered in Subsection 13. 3t 
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